Confirmation of the 1 + Meson Exotics in the rjn System 
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The exclusive reaction iv~p — > rfK°n, rj — * tt + tv~-k° at 18 GeV/c has been studied with a partial 
wave analysis on a sample of 23 492 r/Tr°n events from BNL experiment E852. A mass-dependent fit 
is consistent with a resonant hypothesis for the P+ wave, thus providing evidence for a neutral exotic 
meson with J PC = a mass of 1257 ± 20 ± 25 MeV/c 2 , and a width of 354 ± 64 ± 60 MeV/c 2 . 
New interpretations of the meson exotics in neutral nit system observed in E852 and Crystal Barrel 
experiments are discussed. 

PACS numbers: 13.20.Jx; 13.85-t; 14.40.Cs 
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INTRODUCTION. 



Exotic mesons with J PC = , 1 ^,2^ , . . . do not mix with quark-antiquark mesons and thus offer a natural 
testing ground for QCD. Exotic mesons have been discussed [E H, IE HL HI @> 0> II, H G3 for many years but have 
only recently been observed experimentally. The underlying structure of the negatively charged exotic state with 



J PC = 1 ^ observed in this experiment 11, T2| at 1400 MeV decaying into r]Tr~ is not yet understood. 

Two distinguishing characteristics of the tjtt system make it an excellent one to clarify the properties of this exotic 
state. First, C-parity is a good quantum number in the rjir system, unlike the f]ir~ system. Second, the production 
mechanism for the charge exchange reaction ir~p — > ?/7r°n cannot involve the exchange of an isospin 1 = system 
such as the pomeron. 

The Crystal Barrel experiment [l3| confirmed the existence of resonant structure in the i]ir~ system using stopped 
antiprotons in liquid deuterium in the reaction pn — > ir^ir !]. Later this group analyzed data on pp annihilation at rest 
into and presented evidence for an exotic 1 + resonance in the rjir system with M = (1360 ±25) MeV/c 2 

andT= (220 ±90) MeV/c 2 . 

The rj-ir state has been studied in the GAMS experiment [l5[ with ir~p — > r\ix Q n, r\ — ► 2y, 7r° — > 27 at 32, 38 and 
100 GeV/c. The statistics of the 38 GeV/c data was sufficient so that, using the method of Sadovsky [l6|, they were 
able to present evidence for the exotic 7Ti(1400). 

The VES experiment also observed a peak in the P+ wave of the ijir system near 1400 MeV/c 2 [l7|. In their 
most recent publication [l8| , using theoretical arguments the authors state that the peak can be understood without 
requiring an exotic meson. 

An analysis of E852 data with ir~p — » r]Tr a p, rj — > 27 was recently reported [19j . A bump in the P+ wave of the rjir 
system was observed at M(r]ir°) = 1272 MeV/c 2 with a large width (r = 660 MeV/c 2 ) when fitting all the data using 
the method of Ref. [3). For small t', they observe a width of 190 MeV/c 2 . The authors chose not to claim evidence 
for exotic 7Ti(1400) meson production. 
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The present analysis studies ir~p — > rj7r°n, r\ — > 7r + 7r~7r° at 18 GeV/c in E852. The advantage of this mode over the 
all-neutral final state is that the production vertex is defined by charged tracks. This improves the mass resolution 
as well as the ability to require that the interaction took place within the target. 

EXPERIMENTAL SETUP AND DATA SELECTION 

The data for this analysis were obtained at the Alternating Gradient Synchrotron (BNL USA) in 1995. Using an 
18 GeV/c 7r~ beam interacting in a liquid hydrogen target, a total of 750 million triggers were acquired of which 108 
million were of a type designed to enrich the exclusive final state 7r~p — > 7r + 7r~ A'yn. A total of 6 million events of this 
type were fully reconstructed. The data were kinematically fit [2(| to select events consistent with the 7r~7r + 7r 7r n 
hypothesis (with a confidence level of at least 0.01%) yielding about 4 million events. Of those, 85 228 events passed a 
mass cut enhancing r\ mesons, m(7r _ 7r + 7r°) < 0.65 GeV/c 2 , and 74 549 passed a cut to remove events passing through 
a low-efficiency region in the drift chambers. A final kinematic fit selected 23 492 events for the partial wave analysis 
(PWA), which were consistent with the rpx n, r\ — * tt + tt^tt hypothesis at a minimum confidence level of 1%. 

A strong ij meson signal is observed in this final data sample (Fig. [Iji) with a mass of 539.2 ± 0.3 MeV/c 2 and a 
(resolution-dominated) width of 23.7 ±0.2 MeV/c 2 . The filled regions in the figure indicate the signal region and the 
side-band regions used in the analysis. In the r/ signal region, the signal-to-background ratio is about 4 to 1 for all 
r]7T masses and 5 to 1 for m(riir ) > 1.1 GeV/c 2 . The rjir mass spectrum shown in FigQ}) has two clear peaks: the 
ag(980) and the ag(1320). 

The non-?7 background was estimated as a function of rjir mass using the side-band and signal regions. The 
background fraction varies between 24% and 14% going from lower to higher mass in the region 0.78 < m(r]n ) < 
1.74 GeV/c 2 . 

The experimental acceptance was determined using a Monte Carlo event sample generated with isotropic angular 
distributions in the Gottfried-Jackson frame. The detector simulation was based on the E852 detector simulation 
package SAGEN [l^] . The experimental acceptance was incorporated into the PWA by means of Monte Carlo 



normalization integrals 12}. The acceptance as a function of mass and as a function of t' is flat. 



PARTIAL WAVE ANALYSIS 



The partial- wave analysis (PWA) method described in [12j (see also [2l|, |22j) was used to study the spin-parity 
structure of the r)ir° system in this data set. The PWA was carried out using the extended maximum likelihood 
method separately in each mass bin in the mass region between 0.78 and 1.74 GeV in mass bins of 0.04 GeV for 
< \t'\ < 1.0 (GeV/c) 2 using the likelihood function 

n „ 

ln£cx^ lnI(Qi)- / dffy(fi) /(O). (1) 

i 

Here is the predicted angular distribution, r](£l) is the angular acceptance, and the sum is over the event sample. 

The partial waves are parameterized by a set of five numbers: J PC m e , where J is the angular momentum, P and 
C are the parity and the C-parity of the r/ir system, m is the absolute value of the angular momentum projection 
and e is the reflectivity. We use a simplified notation where each partial wave is denoted by a letter indicating the 
rjn system's angular momentum in standard spectroscopic notation, and a subscript which can take the values 0, +, 
or — , for m e = 0~, 1 + , or 1~ respectively. We assume that the contribution from partial waves with m > 1 is small 
and can be neglected [12[ . 

The amplitudes used are the unnatural parity-exchange waves (UNPW) So, Po, P_, Do, D_, and the natural 
parity-exchange waves (NPW) P+, D + . The NPW waves interfere between themselves as do the UNPW waves but 
the NPW waves do not interfere with the UNPW waves. The P+ wave would be an exotic J PC — 1 h (denoted by 
ir®) if the wave is resonant. 

For each partial wave the complex production amplitudes were determined from an extended maximum likelihood fit 
[22| . The spin 1 /2 nature of the target proton leads to spin- flip and spin-non-flip amplitudes and thus to a production 
spin-density matrix with maximal rank two. The PWA fit p resented in this paper was carried out with the assumption 
that a spin-density matrix of rank one was sufficient [12j. An isotropic incoherent background was included. The 
magnitude of the background was fixed as determined from the side bands. We investigate the quality of the fits by 



comparing the moments of the decay angular distributions H(LM),L < 4 [12j, |22j, of the data with those predicted 
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FIG. 1: (a) Fit of the -K + -K~n° mass distribution in the r/ mass region. There are two entries per event: one for each way 
to assign a 7r° to the r\ decay. The r\ signal region and the side-band regions are shown shaded, (b) The uncorrected t)-k° 
effective-mass distribution for events consistent with the reaction tt~p — » r)ir°n 



by Monte Carlo events generated with the fit amplitudes. We also directly compare the angular distributions for 
cos(6>gj) and tp TJ between the data and those Monte Carlo events. The quality of the fits is good. 

Since natural-parity exchange (NPE) and unnatural-parity exchange (UNPE) amplitudes have different \t'\ depen- 
dences, a fit to the \t'\ distribution using a function of the form N(t') — rii\t'\e~ bl ' t ' + z^e" 62 '* ' was carried out to 
determine the relative contributions of the two exchanges. The fitted parameters are b\ — (7.41±0.08)(GeV/c) 2 , 62 = 
(2.68 ± 0.07)(GeV/c) 2 , and thus the ratio of UNPE and NPE contributions is equal to 0.71 ± 0.03. A value of about 
70% for the ratio of UNPE to NPE at 18 GeV/c is expected based on interpolation of experimental 7r~p charge 
exchange data between 3 and 40 GeV/c in the Regge model [lit ]. 

There is mathematical ambiguity in the description of a system of two pseudo-scalar mesons For our set of 
amplitudes there are eight ambiguous solutions, each of which leads to identical angular distributions. These solutions 
were found analytically starting from one solution by means of the Barrelet zeros method (22j . The eight solutions 
in each mass bin are shown in Fig. [2] as a point for every ambiguous solution. (In some cases the solutions are too 
close together to be visible as separate solutions. This is particularly true in the mass bins at M = 1.56 GeV/c 2 and 
at M = 1.72 GeV/c 2 .) The spread between the various ambiguous solutions dominates the systematic uncertainty in 
the resonance parameter determination (see below). 



MASS DEPENDENT FITS 



To study the resonant structure in the partial waves, we used three different procedures. The first method utilizes 
a Mass Dependent Fit (MDF) of the av erag e solutions in the NPW sector. The PWA results in each mass bin were 



averaged between ambiguous solutions [12[. The mass dependence of the P+ and D + intensities as well as their 



relative phase difference were then fit by relativistic Breit-Wigner (BW) functions (in both the P+ and D + waves) 
with mass- dependent widths and Blatt-Weisskopf barrier factors [l^- The mass and width of the a% are well known 
and were fixed using values^: M = 1320 MeV/c 2 , and V = 120 MeV/c 2 . The width of the a 2 (1320) includes 
the experimental mass resolution. There are three free parameters in the fit of the D + intensity, |-D+| 2 : one for the 
magnitude and two parameterizing the smooth background for the D + wave, as was done in |l2j. 

In the MDF of the |-P+| 2 distribution and the relative phase A$(Z3 + — P + ) there are four free parameters: three 
from the BW function and one for the production phase (assumed constant). The fit was carried out in the mass 
interval 1.1 — 1.74 GeV/c 2 . The resonant hypothesis for D + and P+ waves with a mass-independent production phase 
gives a x 2 /DoF=1.14 for 28 degrees of freedom. The non-resonant hypothesis (no phase variation for the P+ wave) 
gives x 2 /DoF=3.02. It is clear from Fig. [2^ that a single resonant phase for the 02(1320) (dotted line) with a constant 
(non-resonant) P + wave is not satisfactory. 

The contrast between the D + — P + phase variation and the D + phase variation with mass (Fig. solid line) 
shows clearly that the observed phase variation is consistent with interference between two resonant waves. It's worth 
pointing out that the mass dependence of the D + — P + relative phase would be nearly flat if each wave were resonant 
with very similar masses and widths. 

The P + resonant parameters from the fit with the average solutions and the average error matrix [l2l | are: M — 
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FIG. 2: The Partial Wave Analysis (PWA) and Mass-Dependent Fit (MDF) results. The points shown in each mass bin are 
the eight ambiguous PWA solutions, a) the D+ wave intensity; b) the P+ wave intensity; and c) the relative phase between 
the P+ and D+ waves. The lines show the MDF results (method 2). The average PWA solution in each mass bin is plotted 
using grey (yellow) points. The dotted line in (c) is the phase difference if the P+ phase is constant. 



1265 ± 21 MeV/c 2 and T = 411 ± 64 MeV/c 2 . 

The second method was similar to the first except that instead of fitting the average solutions, a large number 
(~ 10 3 ) of randomly chosen combinations of ambiguous solutions in each mass bin were used as input to the mass- 
dependent fit. The obtained distributions of the mass and width of the P + resonance for those fits with acceptable 
values of \ 2 (that is, those with x 2 /DoF < 2) were then fitted by a Gaussian. The mean values of these distributions 
are: M = 1257 ± 20 MeV/c 2 and F = 354 ± 64 MeV/c 2 . The curves shown in Fig. [5] are drawn using these mean 
values. 

The RMS values of these distributions are ctm — 25 MeV/c 2 and or = 58 MeV/c 2 . The systematic errors are 
obtained using these RMS values. Our analysis shows that this spread due to the different ambiguous solutions 
dominates the systematic error. 

The third method used was to carry out a Mass-Dependent Partial Wave Analysis (MDPWA) [T3|. In this procedure, 
an extended maximum likelihood function is generalized to include not only the angular distribution, but also the 
r]ir° mass distribution for each wave. This analysis is free from the problem of ambiguous solutions but it is necessary 
to parameterize the mass dependence of every partial wave (including the UNPWs) and all relative phases. We use 
the same parametrization for the D + and P + waves as in the first two methods. The mass dependence of the UNPW 
waves were chosen to be polynomials of second order with constant phases except for the So wave. The So wave 
was fitted with a BW function using the ao(980) resonance parameters. The MDPWA results for the P+ wave are 
M = 1256 ± 10 MeV/c 2 and V = 319 ± 34 MeV/c 2 , consistent with the results of methods 1 and 2. 

In [l2T | it was shown that a pure D + wave can artificially induce a P + wave due to the effects of finite acceptance 
and resolution. This "leakage" leads to a P+ wave that mimics the D + intensity and phase. In our case, the P+ 
intensity would therefore have an intensity with the shape of the 02 (1320) and a D + — P + phase difference which 
doesn't depend on mass. These features allowed us to include in the MDPWA fit a leakage term with these features. 
We observed that the leakage contribution to the P+ wave from the D + wave is negligible. 

Evidence for a resonance interpretation for the P + wave is primarily the behavior of the D + — P + relative phase 
(Fig. [5b). Since the D + phase variation is well known because of the ai (1320) production, it is clear that the P + 
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phase cannot be constant (see dotted curve in Fig. [2};) and it is well-described by a BW phase variation. 

The ratio of the P + and D + intensities in the range 1.24 < M(t?7T ) < 1.34 GeV is equal to \P + \ 2 /\D + \ 2 = 0.43±0.10. 
This ratio is larger than that for the r)ir~ system, as reported in Ref. [l2| • 

The mass of the neutral exotic 1 h state, decaying into my ) observed here (M = 1257 ± 20 ± 25 MeV/c 2 ) is lower 
than the mass observed in the Crystal Barrel experiment [14[ (M = 1360 ± 25 MeV/c 2 ) by about 100 MeV although 
the results are barely consistent within errors. The width measured here (T = 354 ±64 ±58 MeV/c 2 ) is also consistent 
with that from the Crystal Barrel measurement (T = 220 ± 90 MeV/c 2 ). 

It should also be noted that our result is similar to those obtained from the \ow-t' fits in Ref. [3]. They obtained 
M = 1301±14 MeV/c 2 and T = 190±32 MeV/c 2 in one analysis, and M = 1386±32 MeV/c 2 and T = 363±81 MeV/c 2 
in a fit to the experimental moments. No systematic errors were given for either method. 

The lower mass found in our analysis may be a consequence of interference between the resonant state and back- 
ground in the rjir system, some of which may be from rescattering between the r\ and the ir°. It is also possible that 
two or more resonant 1 ^ states may be present in the r]ir° decay channel in the mass interval mass between 1200 
and 1400 MeV as might be expected if the exotic state is a four quark state. 

The mass for the it® is some 100 MeV below the negatively charged ir^ measured by Crystal Barrel [13J and BNL- 



E852 [ll|,[12|]. Of course, this does not constitute compelling evidence for more than one exotic state in this mass 
region. However, if we assume that there are two ijit states at 1280 and at 1380 MeV respectively, then one expects 
the lower-mass state to have a substantial branching ratio into pn (the reaction tt~ + p — » r/ir + n is mediated by 
p exchange), whereas the higher-mass state would couple to the /2(1270)-7r channel or perhaps to the Pomeron + tt 
channel (the reaction 7r _ + p — > rjn~ + p is produced in part by the exchange of the Pomeron). 



CONCLUSIONS 



Mass dependent fits of the D + and P+ amplitudes and their relative phase using three different methods (described 
above) all lead to the conclusion that the P + wave is well-described by a resonance hypothesis and is inconsistent with 
having a constant phase. The resonance parameters for the observed 7r° are given by M = 1257 ± 20 ± 25 MeV/c 2 
and r = 354 ± 64 ± 58 MeV/c 2 . Here the first error is statistical and the second is systematic. (We have chosen to 
take the resonant parameters and errors from method 2.) This result, together with the previous results from Crystal 
Barrel 14 1 and E852 [l|| provide strong indications for one or more spin-exotic mesons near 1400 MeV/c 2 decaying 
to r]ir°. 
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